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STM measurement of single spin relaxation time in superconductors. 
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Localized spin states in conventional superconductors at low temperatures are expected to have 
long decoherence time due to the strong suppression of spin relaxation channels. We propose a 
scanning tunneling microscopy (STM) experiment allowing the direct measurement of the decoher- 
ence time of a single spin in a conventional superconductor. The experimental setup can be readily 
applied to general-purpose spin-polarized STM and to local spin relaxation spectroscopy. A possible 
extension of the setup to a quantum information processing scheme is discussed. 



When a magnetic atom is placed in a conventional s- 
wave superconductor, the exchange interaction between 
the localized spin and the conduction electrons leads 
to formation of spin polarized scattering resonances in- 
side the superconducting energy gap These reso- 
nances manifest themselves as sharp peaks in the den- 
sity of states localized near the impurity. The spin- 
relaxation time for such states is governed by the cou- 
pling to the environment. At low temperatures, this 
coupling is strongly suppressed due to the gapped na- 
ture of quasiparticles in s-wave superconductors. The 
weak spin-relaxation channels related to magnetic dipole- 
dipole interaction with the nuclear magnetic moments 
can be further reduced by choosing a superconductor 
with zero nuclear spin. It is therefore expected that 
long spin-relaxation times are achievable when a local 
moment is placed in s-wave superconductor. In this Let- 
ter we propose an experimental setup capable of deter- 
mining the single-spin relaxation time for such defects. 
Provided that the measured spin-relaxation time is suf- 
ficiently long, superconducting tips ending with a single 
magnetic atom can be used for low-temperature spin- 
polarized STM. We also discuss a possible application 
of the spin-polarized local states in superconductors in 
quantum information processing schemes. 

The tunneling current between two spin-polarized 
states depends strongly on the angle between the respec- 
tive spin polarizations. The time-dependence of the tun- 
neling current can hence be used to extract the dynamics 
of the individual spins, assuming they are weakly cou- 
pled. This is the basis for our method of determination 
of the single-spin relaxation time in superconductors. In 
the proposed experiment, both sides of the tunnel junc- 
tion are identical superconductors, each containing one 
magnetic atom. The spin-polarized signal will be opti- 
mized when the impurities on the two sides of the junc- 
tion are spatially aligned; this is most easily achieved 
in the STM setup, when one of the magnetic atoms is 
placed directly on the tip. Superconducting (for exam- 
ple, Nb) STM tips have already been successfully used 
H and the STM capabilities for the single-atom manip- 
ulation H make the fabrication of the tip ending with a 
single magnetic impurity atom feasible. 



To theoretically predict the dependence of the I-V 
characteristic on the angle between the impurity spins on 
the tip and on the surface, we use the T-matrix approach, 
combined with the tunneling Hamiltonian formalism. 
We model the impurity by introducing a local pertur- 
bation to the Bardeen-Cooper-Schrieffer (BCS) Hamil- 
tonian describing a homogeneous superconductor. In 
the Nambu formalism || , the impurity Hamiltonian can 
be written in terms of two-component Nambu spinors, 
*k = [c kT cL k ,] T , as @ 
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Here V — Wf + Uf 3 , f is the 2x2 unit matrix, 
Ti(i — 1,2,3) are the Pauli matrices. The impurity 
strength is characterized by the exchange interaction 
with the conduction electrons, W, and the on-site poten- 
tial, U, or the corresponding dimensionless parameters 
w = ttNqW and u — itNqU ', where No is the density of 
states at Fermi level in the normal state. We can then 
obtain the expression for the complete finite-temperature 
Green's function (GF) in the presence of single impurity, 

+ g^{v-iu J )f(iu J )g^{-v'-iu J ), (2) 

in terms of the GF for the homogeneous SC g(°\r;iu>) 
and the T-matrix f(iw) = [V^ 1 - £ (0) (r = 0; We 
approximate the finite-temperature GFs for the homoge- 
neous SC in the real-space representation by 
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where 2 A is the SC energy gap and kp is the Fermi 
momentum. By using expression ([|) we neglect the 
exponential decay of the GF for the homogeneous su- 
perconductor with r = |r| on the length scale of the 
superconducting coherence length (typically ~ 1000 A 
for conventional SCs) and the effects of the finite band- 
width. Both approximations are justified, because we 
are only interested in short-length (of the order of few 
Angstroms) and low-frequency (of the order of SC gap) 
behavior of the GF. We also neglect the suppression 
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of the SC order parameter in the vicinity of the impu- 
rity, because this approximation does not significantly 
affect the spectral functions 0. Using this set of ap- 
proximations, we calculate the GF for a superconduc- 
tor with an impurity and found the expressions for the 
spectral functions, corresponding to its diagonal elements 
Aa(r, r'; e) = — (l/7r)ImG^(r, r'; ±e). Here the plus sign 
in the argument of the retarded GF G^(r, r'; e), obtained 
by analytic continuation of (0), corresponds to i = 1 
and minus sign to i = 2. Spectral functions An have 
a particle-hole-symmetric continuous part at frequencies 
|e| > A and a S- function singularities on the symmetric 
energies, 



to = (-)' 



vA 



(4) 



where a — (1 + u 2 — w 2 )/ (2w). 

In the tunneling Hamiltonian formalism, the quasipar- 
ticle tunneling current can be expressed as a convolution 
of the spectral functions on different sides of the junc- 
tion JtJ. In our case, the impurities on the tip and on 
the surface (we will use sub- or superscript "1" to refer 
to the tip properties and "2" for the surface properties) 
may have different spin orientations. To account for this, 
we introduce angle 9 between the spin directions and ex- 
press the fermion creation and annihilation operators on 
the surface in the spin basis with the quantization axis 
along the tip impurity spin. The result is 
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Here f, for example, stands for the "up" spin orientation 
in the surface spin basis and f is the "up" spin in the tip 
spin basis. Then, we come to the expression for the total 
quasiparticle current 



/(F) = (/n + / 22 ) cos 2 
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7 21 )sin 2 -, (6) 



with the partial currents, Iij(V), given by 
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x A^{v i y i -e)Af>{r 2 y 2 -e-eV)de. (7) 

Here e is the electron charge, ®{e,eV) — np(e — eV) — 
np(e) is the difference of the Fermi-Dirac distribution 
functions on the different sides of the junction, T riiF2 is 
the tunneling matrix element, and the summation is per- 
formed over all lattice sites on the tip (ri, r[) and on 
the surface (r2, r' 2 ). Expression ^ is for general elec- 
trode geometry. We can further simplify it by assuming 
that the tunneling happens from a single point on the 
tip (the impurity atom). The tunneling matrix element 
then takes the form T{r) — Texp(— \Jr 2 + z 2 /d) that de- 
pends on the lateral distance r from the projection of the 



tip onto the surface, tip-surface distance z, and the char- 
acteristic decay length d, related to the material work 
function (typically about 0.5A). Detailed calculations for 
this and other geometries will be presented elsewhere. 

To calibrate the theoretical parameters for the impu- 
rity and the superconducting host, we consider tunneling 
experiment from a normal STM tip into the surface of 
the SC (Nb) with the magnetic adatoms (Gd or Mn) on 
it 0. We can qualitatively reproduce experimental tun- 
neling spectra by tuning the impurity parameters w and 
u. Results of this calculation are presented in Fig. [j]. In 
what follows, we adopt parameters characteristic of Mn 
adatom in Nb: w = 0.9, u = 1.0, and A = 1.48 meV. 
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FIG. 1. The calculated change in the dl/dV spectrum 
due to the presence of Mn (a) and Gd (b) adatoms. Impurity 
parameters were taken to be w — 0.9, u = 1.0 for Mn and 
w = 0.3, u = -0.3 for Gd. The surface is Nb with A = 1.48 
meV at T — 3.85 K. Tip-surface distance is z = 5 A and 
d = 0.5 A. 



When the impurities are present both on the SC tip 
and on the SC surface, a number of new features appear 
in the I-V characteristic. For example, when voltage first 
becomes large enough (in absolute value) for the intragap 
state of the tip to overlap with the continuous part of the 
spectrum of the surface, the abrupt increase (or decrease) 
in current results. The most remarkable feature of these 
I-V characteristics, however, is their strong dependence 
on the relative orientation of the impurity spins, i.e. on 
the angle 9. The reason for this strong dependence is 
the existence of the contribution to the current, resulting 
from the overlap between the intragap states of the tip 
and of the surface (at certain values of voltage). From the 
expression for the position of the intragap states, Eq. (Q), 
it can be seen that these peaks (generally there are four 
of them) appear at voltages 

«iAi , a 2 A 2 



eV = ±- 



± 



(8) 



For the identical impurities in the identical superconduc- 
tors, (a\ = a 2 ), Ai = A 2 , at 9 = these peaks vanish 
because the corresponding intragap states, contributing 
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to In and 1 22 overlap only at zero voltage and hence there 
is no tunneling current between them. Fig. |^ presents the 
I-V characteristics calculated for various values of angle 
9. The (5-function peaks due to the tunneling between 
the intragap states appear for non-zero 8 at V = ±1.63 
mV and become more intensive with increasing 8. It is 
this strong 8 dependence that allows us to formulate the 
experimental procedure for the measurement of the de- 
coherence time of the intragap impurity spin-polarized 
states in the SC. 
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FIG. 2. Calculated I-V characteristic of the Nb super- 
conducting tunneling contact for different values of angle 8 
between Mn spins on two sides of the junction. The current 
is normalized by Iq = <7oAi/e, where ao is the normal-state 
conductance of the tunneling junction. The curves are offset 
by Jo and the 5-functions at V = ±1.63 mV (dashed vertical 
lines) are broadened for clarity. 



The simplest way to measure the spin decoherence time 
is to monitor time dependence of the tunneling current, 
I(t), for applied voltage fixed at the resonant value (V = 
±1.63 mV in our calculation). The current-current cor- 
relation function, s(t) = ((i(t)i(o))-(i) 2 )/((i 2 )-(i) 2 ), 

can be used to directly extract the spin-relaxation time. 
The theoretical prediction for S(t) is presented in Fig. [|. 
Another way to extract spin dynamics is to first polarize 
the spins by external magnetic field (larger than ksT/ fis 
but smaller than H c i) and then observe the relaxation of 
the current after the field is removed. The effect of the 
magnetic dipole-dipole interaction between the impurity 
spins can be studied by measuring the current relaxation 
at different tip-surface distances. It is possible, however, 
that the time resolution of the experiment is insufficient 
to resolve the times associated with the spin dynamics; in 
this case, the measured current will be a constant equal 
to the average value in Eq. (0). 




FIG. 3. Current-current correlation function 

S{t) = «/(i)/(0)) - </) 2 )/«/ 2 ) - </) 2 ), calculated for un- 
corrected diffusion of spin orientations. Assuming that each 
spin diffuses over the unit sphere with the diffusion coefficient 
A, the relative angle between the spins evolves according to 
stochastic PDE, d9 = A cot(0)di/2+v / Ad£, where £ is the nor- 
mal random variable with zero mean and variance (d£ 2 ) = dt. 
This correlation function is well approximated by the usual 
expression for diffusion, exp(— At). 



Generally, the same technique may be used to measure 
the relaxation time of any local spin state on a conduct- 
ing substrate, as long as this relaxation time is shorter 
than the relaxation time of the impurity-induced intra- 
gap state in the SC tip, e.g. Kondo spin dynamics on a 
metallic surface 

If the relaxation time of the intragap spin-polarized im- 
purity state in a SC proves to be sufficiently long, a SC 
tip ending with a magnetic impurity can be used to ob- 
tain the spin-contrast STM images. Currently, to achieve 
the spin contrast in the STM measurements one has to 
use ferromagnetic tips PJTc|]. These STM experiments 
have already produced a number of interesting results, 
like the direct imaging of the two-dimensional antiferro- 
magnetism E|. Unfortunately, in such setup, the increase 
of the tunneling area of the tip due to the restrictions on 
the tip material leads to the degradation of the spatial 
resolution, compared to conventional STM. In addition, 
strong magnetic fields produced by the tips may affect 
the measurements by inducing polarization of the local 
spin structure, or in the case of superconducting sam- 
ples, cause a local SC-to-normal transition. Our pro- 
posed setup would be the next step in the spin-polarized 
STM, which would allow the spin-contrast imaging at low 
biases with an atomically sharp STM tip. In this setup, 
the effect of the magnetic field due to the magnetization 
of the tip is minimized since a single impurity spin would 
be sufficient to achieve the spin contrast. As an example, 
in Fig. [|, we present a theoretically calculated map of 
the spin-contrast signal, 
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in the vicinity of surface impurity. Currents in Eq. (|9|) 
correspond to the parallel (Ijf ) and antiparallel (Jf j.) ori- 
entations of the impurity spins on the tip and on the 
surface. 
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FIG. 4. The topographic map of the log of the 
spin-contrast STM signal P in the vicinity of the magnetic 
impurity. Rapidly decaying Friedel oscillations are clearly 
seen at large distances. For this calculation we have taken 
the same values of parameters as for Fig. ^ and fixed the bias 
voltage at V = 1.63 mV. 

Another application of the spin polarized states near 
magnetic atoms in superconductors is related to quan- 
tum information processing. It is motivated by the ex- 
pected long relaxation times of spins in s-wave super- 
conductors and the ease of read-out of the spin states 
using the spin-polarized tunneling techniques discussed 
above. At temperatures much lower than the SC tran- 
sition temperature, the quasiparticle contribution to the 
spin relaxation will effectively disappear. The dipolar re- 
laxation due to interaction with the lattice nuclear spins 
can be also significantly reduced by using SC with non- 
magnetic nuclei, e.g. non-magnetic isotopes of Pb or Sn. 
A local spin in the s-wave superconductor is isolated from 
environment to the same extent as a local spin in a semi- 
conductor with the band gap equal to 2A. An advantage 
compared to current proposals based of semiconductors, 
e.g. phosphorus in silicon [jy], is that the substrate is 



conductive and hence one can do local tunneling exper- 
iments to interrogate local spins that store and process 
quantum information. The details of the proposed quan- 
tum information procession architecture based on the lo- 
cal spins in superconductors will be discussed elsewhere 

& 

In summary, we have proposed an STM experiment, 
which will allow the measurement of the decoherence 
times of the single spins in superconductors. For suffi- 
ciently long measured spin-decoherence time, the setup 
can be applied to the spin-resolved STM measurements. 
Also, the possibly long decoherence time of the single- 
impurity spin states in superconductors, combined with 
the ease of readout and manipulation, make such states 
a favorable candidate for a qubit in quantum computing 
applications. 
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